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Abstract Secular and seasonal variations of the coefficient J2 of the geopo- 
tential are studied from the analysis of laser measurements of distances to the 
geodetic satellites LAGEOS 1 (1988-2003) and LAGEOS 2 (1992-2003). It 
n , ' is confirmed that beside the well-known annual variations with the amplitude 

[tJ ' K, 2.5 X 10"^" there also exist very significant semi-annual variations of a com- 

parable amplitude. Phases of these two modes are such that the total effect 
may be described as a sharp postive splash of J2 in August and considerably 
smaller variations in the rest part of year. The adopted theoretical value of the 
O ' so-called dynamical Love number fcj (a scale factor of near-diurnal oscillations 

aroused by the Earth's fluid core in the coefficients c\ , s^ of the geopotential) is 
jrt ' improved applying a simple close form of these oscillations expressed in terms 

of the differential angular velocities vi, V2 of the fluid core. It is shown that 

this form is equivalent to the standard one of Fourier series in which such os- 

^ , dilations usually are referred to as a frequency-dependence of the Love number 

1/^ ' fc2. The derived estimate /cf — 0.0595 ± 0.0007 statistically differs from the 

^^ , theoretical value fcf — 0.063. Out-phase oscillation of fcj with the period about 

^^ \ 18.6 years and the amplitude 0.0064 ± 0.0011 is detected giving evidence of 

large dissipation in the fluid core. The estimated secular trend J2 (commonly 



C"*~- ' interpreted as the effect of the so-called Post-Glacial Rebound) appears twice 



less than the value recommended by the standards of the International Earth 
Rotation Service (lERS) but agrees with last findings of other authors. 

Keywords: Earth's satellites, geopotential, tides, Love numbers, SLR ob- 
servations 



1 Introduction 

Starting point of this research was studying the near-diurnal oscillations dc\, 
ds\ caused by the Earth's fluid core in the coefficients c\, s\ of geopotential. 
This effect (interpreted as the frequency-dependence of the Love number fc2) is 
presented by lERS standards in the form of the sum of trigonometric harmonics 
that depend on the fundamental arguments of the nutation theory (McGarthy, 
Petit, 2004). It may be shown that the amplitudes of these harmonics are 
proportional to the so-called dynamical Love number fc^. In Appendix, the 



expression for the corrections dc2, ds\ is derived in the simple close form: 



d4\ ^_^-ij M f ^1 



J^l^ , > (1) 



where vi , V2 are the equatorial components of the angular velocity of the fluid 
core relatively to the Earth as a whole, J2 is the coefficient of the second zonal 
harmonics of the geopotential, lo is the Earth angular velocity, kg — 0.93831 is 
the so-called secular Love number. It will be also shown that this presentation 
is equivalent to the standard form of trigonometric series for the frequency- 
dependence of the Love number k2 ■ 

The parameter fc2 plays important role in two geodynamical problems. Firstly, 
it is rotation of the deformable Earth in the framework of the standard SOS 
model (Moritz, Muller, 1987), and secondly, dynamics of the Earth's satellites. 
In short, perturbations of rotation of the deformable Earth depend on the so- 
called compliance £, that describes the tidal response of the moments of inertia of 
the Earth on the differential rotation of its fluid core. This parameter may be ex- 
pressed in terms of the dynamical Love number k2 by the relation ^ = efc2 /fc^ in 
which e is the dynamical flattening of the Earth. Unfortunately, while process- 
ing Celestial Pole positions provided by geodetic VLBI observations it proved 
impossible to separate A;| from the correction to the ratio of the moment of in- 
ertia of the fluid core to that of the Earth as a whole (Shirai, Fukushima, 2001), 
(Krasinsky, Vasilyev, 2006). On the other hand, the near-diurnal variations of 
the coefficients c\ , s\ of the geopotential bring about accumulating secular ef- 
fects in the node of satellite orbits and thus one may hope to estimate ^2 from 
satellite laser ranging (SLR) observations of geodetic satellites. As the adopted 
value of fc2 is purely theoretical, confronting it with experimental data is quite 
an actual problem. The value derived in this paper from the analysis of SLR 
data of the geodetic satellites LAGEOS 1 and LAGEOS 2 is, to our knowledge, 
a first experimental estimation of this parameter. 

^From the same analysis of the SLR data, we determine also corrections to 
the coefficient J2 of the geopotential for studying the time-variability of this 
coefficient. In such a study, obtaining the secular trend J2 is a problem of 
primary importance as the lERS standards recommend to take this effect into 
account whenever numerical integration of the equations of artificial satellite 
motion has to be done. At present, two centers of analysis (University of Texas 
and Goddard Space Flight Center) are monitoring the geopotential on a regular 
basis. Recent publications of the first group (Cheng, Tapley, 2008) and the 
second one (Lemoine et al., 2006) have revealed serious discrepancies between 
the derived values of the secular trend: J2 = —2.6 x 10^^^/ycar and J2 = —1.3 x 
10"^"'^ /year, correspondingly. So, the independent study of J2 documented in 
Section 3.3 seems to be actual. 



2 Observations and their processing 

We have used the laser observations of the geodetic sateUites LAGEOS 1 (1988- 
2003), and LAGEOS 2 (1992-2003) taken from the server' ftp://cddis.nasa.gov\ 
The observations of each year were combined into one series and treated simuha- 
neously. In Tab.l, the following information on the observations is given for each 
year: their total number Nobs, the total number Ndei of rejected observations, 
the Weighted Root Mean Square (WRMS) error of one-way distances to the 
satellite (in millimeters) . Six orbital elements of the satellite under study and 
two dynamical empirical terms (along- and cross-track perturbing accelerations) 
were estimated weekly as local parameters. When processing the observations 
of each year, these local parameters were determined simultaneously with the 
global ones which are the coordinates of observing stations and the value of 
the dynamical Love number fcj. In addition, on each monthly interval (more 
exactly, on the interval of four weeks) a correction to the coefficient J2 of the 
geopotential also was determined for further study of its time variations. It 
is known that SLR observations made in different stations are of quite various 
quality. Properly down-weighted or rejected, the observations of poor quality 
do not distort the resulting estimates. In fact, only few stations of the last 
generation provide the overwhelming volume of high precision data that affect 
the results. 

As an example, from the typical set of the observations of LAGEOS 1 of the 
year 2000, the following parameters have been derived: 



• dynamical Love number /cj, 

• monthly corrections to J2 (13 unknowns), 

• weekly corrections to six elements of the satellite and two empirical pa- 
rameters (52 X 8 = 416 unknowns), 

• yearly corrections to the coordinates of 99 stations (295 unknowns). 

Thus, the total number of the estimated parameters for this year is 626. 
Values of the longitudes and latitudes of several fiducial stations (which provide 
the bulk of accurate observations) have been fixed to prevent degeneration of 
the normal matrix due to correlations with the elements of orientation of the 
satellite orbit. The software complex ERA for ephemeris and dynamical as- 
tronomy was applied while constructing the numerical theories of the orbital 
motion, calculating the theoretical distances to the satellites and processing the 
condition equations. The current DOS and Windows versions are available from 
the anonymous FTP server quasar. ipa.nw.ru/incoming /era. All calculations are 
carried out in accordance with the recommendations of lERS 2003 conventions 
(McCarthy, Petit, 2004). The only difference in the dynamical equations of 
the satellite motion is making use of the close form ([T|) for the corrections dc\^ 
ds\ caused by the fluid core. We will show that this approach is equivalent 
to that recommended by lERS but is more convenient because the functional 
dependence of the perturbing forces on fc2 is presented in explicit form. 



Table 1: WRMS errors a of the post-fit residuals (nun) 







LAGEOS 1 






LAGEOS 2 




Year 


Nohs 


Ndel 


a{mm) 


Nobs 


Ndel 


WRMS (mm) 


1988 


67414 


3058 


22 








1989 


64030 


3411 


22 








1990 


77951 


2721 


22 








1991 


54981 


3518 


22 








1992 


59715 


6073 


24 


13683 


863 


21 


1993 


81398 


7665 


22 


78873 


7128 


22 


1994 


65370 


4990 


27 


66827 


4557 


21 


1995 


53463 


4353 


24 


52153 


3183 


20 


1996 


52725 


4196 


23 


52244 


7904 


23 


1997 


50404 


3336 


22 


53516 


3688 


17 


1998 


65010 


5361 


19 


60819 


3298 


17 


1999 


71788 


4693 


23 


66520 


3415 


18 


2000 


61888 


4122 


23 


63158 


2362 


19 


2001 


71839 


4060 


19 


70449 


3083 


19 


2002 


70780 


2503 


21 


64011 


2217 


18 


2003 


74969 


2379 


19 


74613 


2798 


19 



3 Results of data processing 



3.1 Dynamical Love number k: 



Applying the method described above, we have obtained the set of 16 yearly es- 
timates of fcf from the observations of LAGEOS 1 (1988-2003) and 12 estimates 
from LAGEOS 2 (1992-2003). They are presented in Fig.l (the black circles 
correspond to LAGEOS 1, the hollow ones to LAGEOS 2). Generally, all the 
estimates are in a very good accordance (the only exception is the LAGEOS 1 
result for the year 1995 which evidently drops out and thus has not been used). 
From this set, the following averaged value of the dynamical Love number is 
derived 



A;2 = 0.0595 ± 0.0007 

that statistically differs from the theoretical a priory value k!^ = 0.063 given in 
(Moritz, MuUer, 1987). 

Fig.l demonstrates also an out-phase oscillation with the period 18.6 years 
and the amplitude 0.0064±0.0011. Probably, that is evidence of large dissipation 
in the fluid core. Note that parameters of the dissipative effects in the fluid core 
recommended by lERS standards are purely theoretical, and their numerical 
values have not yet been verified by observations. 



Figure 1 : Yearly estimates of dynamical Love number k^ 



o 


075- 
050- 


* 










* 










# 


^ 


♦ 


* 


* 


o 




* 


* 


4- 




o 


* 
■o- 


o 


f? 


4- 

o 






























H» 
















o 


035- 
































n 


ooo_ 

































































88 89 90 91 92 



For checking a possible dependence of fc2 on its starting value, the observa- 
tions of an arbitrarily chosen year were processed with the zero starting value 
of A;2 . After several iterations, the resulted estimate appeared the same. 

Note that in the methodically correct approach, the standard "static" Love 
number ^2 should also be estimated from the analysis of SLR data. At present, 
the most accurate and reliable estimate of ^2 raight be obtained from the anal- 
ysis of geodetic VLBI observations. Using these data, the following consistent 
estimates have been derived: k2 — 0.2788 ± 0.0011 (Shirai, Fukushima, 2001) 
and fe = 0.27272 ± 0.00036 (Krasinsky, Vasilyev, 2006) which are significantly 
less than the adopted value ^2 = 0.3. The error of this value may distort results 
of processing the SLR data. 



3.2 Periodic variations of the coefficient J2 of geopotential 

Fig. 2 presents the monthly corrections to J2. The noticeable positive trends 
means that the a priory value J2 = —26 x lO^^^/year used in our calculations 
should be considerably diminished. When comparing Fig. 2 with the analogous 
plots in (Cheng, Tapley, 2004) and (Cox et al, 2005) it is necessary to bear in 
mind that Fig. 2 presents corrections to the assumed lERS value of J2 while those 
in the cited works correspond to the zero value of it. More details concerning 
the value J2 are discussed in the next section. 

Methodically correct estimates of J2 may be derived only simultaneously 
with determination of cosine and sine coefficients of the main harmonics of the 
trigonometric polynomial approximation of J2 . These are the annual harmonics 



Acos, ^sin, the semi-annual harmonics A, 



(0.5) 



4(0.5) 



the harmonics A^^^, A^„ 



connected with the period 18.6 year of revolution of the lunar node J7, and the 
empirical harmonics ^cos, ^gin with the period of 3 year. We also estimated 
the correction dJ2{0) to the value J2 for 2000.0. Four solutions are obtained 
for different sets of the estimated parameters (see Tab. 2). The first column of 
Tab. 2 specifies parameters under estimation, others give their estimated values. 
For the secular trend, there are given both the correction dJ2 to the a priory 
value —26 x 10^^^/year, and the full value J2. 

Tab. 2 demonstrates that the estimates of the annual and semi-annual har- 
monics are robust and only slightly depend on the solutions. For the annual 
oscillations, the results can be compared with those by other authors. The full 
amplitude of the annual oscillations is (229 ±23) x 10~^^ in accordance with the 
value 300 x 10^^^ taken approximately from Fig.l of (Cox et al., 2005). Similar 
results are obtained in a number of other works; for instance, in (Cheng, Tapley, 
2004) the seasonal amplitude is 290 x 10~^^ (no error is given). The phases of 
the oscillations also agree: the negative minimum in January and the positive 
maximum in July. The amplitudes of the semi-annual oscillations are rather sig- 
nificant and comparable with those of the annual ones. The derived amplitudes 
of the semi-annual variations are consistent with the results published by other 
authors. For instance. Table 1 of paper (Nerem et al, 2006.) gives the estimated 
cosine and sine amplitudes of the semi-annual variations as —2.41 x lO""'^" and 
0.78 X 10"-'^° in accordance with our findings (see Table 2). 

The solid curve in Fig. 2 presents Solution 4 (with maximal number of esti- 
mated parameters). Fig. 2 shows that the seasonal variations cannot be described 
as a simple harmonic oscillation but have more complex structure. Namely, they 
keep near constant value from the beginning of every year till August when a 
sharp maximum takes place, and then decrease up to the end of the year. The 
plot given by Fig. 2 is very similar to the intra-annual variations of J2 presented 
by Fig 3(a) of (Moor et al. 2006) and derived from Lageos 1,2 data for the 
time span 1998-2004. Such behavior is a consequence of the large value of the 
semi-annual amplitude. This feature of the intra-annual variations qualitatively 
agrees with Fig. 3 of paper (Lemoine et al. 2006) where they are presented for 
the interval 2000-2007 being derived from SLR data for 9 satellites combined 
with DORIS data. Note that Fig. 3 in (Lemoine et al. 2006) presents variations 
of the coefficient C20 = —^2- It is interesting that Fig. 2 of the same paper does 
not demonstrate the fine structure of the intra-annual variations, probably be- 
cause time-intervals between the consequent samples are two months (though 
the caption indicates one month, probably erroneously). With so low time res- 
olution, the fine structure of the intra-annual variations cannot be discerned. 
Supposedly, this fine structure is due to the seasonal asymmetry of ice cover de- 
veloping in the southern and northern semi-spheres. For a more detailed study 
of these variations, Fig. 3 presents them vs. JD (mod 365.25) after removing the 
secular and long-periodic variations taken from Tab. 2. Beside the maximum in 
August, this plot demonstrates also and 1 ong-periodic variations taken from 
Tab. 2. Beside the maximum in August, this plot demonstrates also two small 
minima at April and November. Here only the LAGEOS 2 results have been 



Figure 2: Corrections to the nominal J2 (for the adopted J2) in 10 ^° 




used as the more accurate ones. The bars present statistical errors obtained 
after averaging 13 groups of the monthly estimates J2 for all the years. 

In Solutions 3 and 4, the sine and cosine components of the oscillation with 
the period 18.6 years have also been estimated (with the resulted WRMS errors 
considerably decreased). The nature of the empirical oscillation of the 3-year 
period estimated in Solution 4 is unclear. We tried several empirical long- 
periodic harmonics but statistically significant amplitudes were found only for 
the 3-year period. 

Fig. 2 demonstrates that the positive outliers correspond to the August spikes. 
It means that the values of these sharp spikes vary with time and cannot be 
presented by the simple model as the sum of the annual and semi-annual oscil- 
lations. Probably, this effect should admit some geophysical interpretation. 

In the last two lines of Tab. 2, the WRMS errors of the residuals are pre- 
sented (both for LAGEOS 1 and LAGEOS 2). Note that the WRMS errors for 
LAGEOS 2 are considerably smaller than those for LAGEOS 1. Formal erors 
of J2 appeared to be by order less than the range of the fluctuations of (iJ2, and 
so could not be presented by error bars in Fig. 2. 



3.3 Secular trend J2 and the post-glacial rebound 

True value of the trend J2 is quite important as for practical applications as for a 
geophysical interpretation. Indeed, the estimate J2 = —2.6 x 10^^^/year is still 
recommended by lERS standards (McCarthy, Petit, 2004) as experimentally 
confirmed effect of the post-glacial rebound, despite that this trend has greatly 
diminished or even completely vanished for the SLR data after 1998 (Cheng, 



Figure 3: Intra-annual variations of J2 vs. JD (mod 365.25) (in 10 ^^) 
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Table 2: Models of variations of J2, four solutions (in 10 ) 





1 


2 


3 


4 


dJ2{0) 


218(16) 


215(14) 


103(31) 


77(33) 


rfJ2/year 


33(3) 


35(3) 


18(7) 


13(7) 


J2/year 


7(3) 


9(3) 


-8(7) 


-13(7) 


^cos 


-172(19) 


-174(17) 


-184(15) 


-183(15) 


^sin 


-151(19) 


-150(16) 


-164(15) 


-168(15) 


4(0-5) 

^T-COS 




66(16) 


67(15) 


66(15) 


4(0-5) 




167(17) 


168(15) 


168(15) 


^cos 






154(39) 


181(40) 


^sin 






-102(20) 


-98(20) 


4(3) 








-43(15) 


.(3) 








-39(16) 


WRMS, LAGEOS 1 


285 


255 


245 


235 


WRMS, LAGEOS 2 


195 


155 


135 


135 



Tapley, 2004). For instance, Fig. 3 in (Chao, 2003) shows that the accumulated 
error of J2 becomes intolerably large at present epoch if one applies the adopted 
value of J2- In recent publication (Cheng, Tapley, 2008), the authors insist 
that this value is still valid and explain the clearly seen deviations of J2 from 
observations after 1998 by some decadal variations, making stress on the effect 
of El Nino. On the other hand, the independent analysis carried out in Goddard 
Space Flight Center (Lemoine et al, 2006) demonstrates that the time-behavior 
of J2 was rather stable since 1976 up to 2007 with no noticeable change of it after 
1998. The resulting value J2 = —1.32 x 10^^^/year is twice less than that of the 
lERS standards based on the works by Cheng and Tapley. Our analysis of SLR 
data confirms this conclusion. In Solution 4 (maximal numbers of the estimated 
harmonics in J2) we have obtained the estimate J2 = (—1.3 ±0.7) x 10^^^/year 
which is consistent with the result of (Lemoine et al, 2006). Unfortunately, 
due to the comparatively short 16-yearly time-interval of the data used, the 
correlations of the amplitudes of the long-term variations with the secular trend 
are rather large. That is why the statistical error of our estimate of J2 is 
comparatively large. In any case, we may state that our value of J2 rules out 
the rate —2.6 x 10^^^/year recommended by lERS standards. 

4 Concluding notes 

1. From the practical point of view, the most important result of this study is 
the conclusion that the derived secular trend J2 is consistent with the result of 
(Lemoine et al, 2006) confirming that the value of this trend recommended by 
lERS standards is twice too large. To our opinion, poor predictability of the 
variations of J2 due to inter-annual oscillations of unclear origin necessitates 
monitoring of the current values of J2 and disseminating the results for practical 
applications, just as at present it is done for the Earth's orientation parameters. 

2. The intra-annual variations of J2 are not a simple harmonic oscillation; 
see Fig. 3. Probably, the sharp spike in August is due to asymmetric seasonal 
developing of the ice regime in the southern and northern semi-spheres. This 
problem deserves thorough study and a geophysical interpretation. 

3. The derived estimate fcg — 0.0595 ± 0.0007 statistically differs from its 
theoretical value fc2 = 0.063. To our knowledge, that is the first experimental 
determination of this parameter. 

4. For independent calculation of various tidal effects caused by the fluid 
core, it is desirable to provide users by numerical values of the angular velocities 
vi , V2 of the fluid core. To avoid somewhat bulky analytical manipulations 
directly in the text of present paper, derivation of the analytical expressions for 
vi , V2 is transferred into Appendix. 



5 Appendix. Dynamical Love numbers 

It is known that the positions of sites on the Earth's surface are displaced by 
the pole tides of two types. Those of the first type are caused by the motion 
of the pole of the Earth as a whole, the corresponding tidal amplitudes being 
proportional to the Love numbers h2, h- The pole tides of the second type 
are caused by the differential rotation of the fluid core relatively to the mantle; 
their amplitudes are proportional to the so-called dynamical Love numbers ft-2: 
12- The pole tides of the both types also contribute to the tesseral coefficients 
C2, si of the geopotential. Putting aside the Chandler's free wobble of the pole, 
these contributions are of near-diurnal periods with the amplitudes proportional 
either to the static " potential" Love number k2 or to the dynamical Love number 
/cg, the latter effect being the largest. A theory of the dynamical Love numbers 
is given in monograph (Moritz and MuUcr, 1987). 

As the fluid core brings about diurnal oscillations of the coefficients c^ and 
S2 of the geopotential, it gives rise to rather significant long-term and secular 
perturbations in the elements of the observed satellites (mainly in the node) thus 
making it possible to improve the adopted theoretical value of the scaling factor 
fcg of these perturbations. It is commonly assumed that the fluid core rotates 
relatively to the mantle with the angular velocity v = (wi,W2, 0), the equatorial 
components vi, V2 being given by the theory of the Earth's rotation. We denote 
oj = (a;i,a;2, W3) the vector of the angular velocity of the Earth as a whole. For 
calculating the potentials dW^ induced by the centrifugal acceleration of the 
differential rotation of the fluid core, the velocity of any point r in these two 
domains within the Earth (the mantle and fluid core) may be presented either 
as w X r or as (cj -I- w) X r, respectively. Then the centrifugal acceleration W 
within the mantle may be presented as 

W — —Zd X (ZU X r) = —uj{r, uj) + ruj'^ , 

while in the fluid core the corresponding expression has the form 

W = -{uJ + v) X [(ZJ + v) xr] = -{uJ + v) [{r,uj) + {r,v)] + r |uJ + iJp 

with the standard notations {p, q) and p x g for the scalar and vectorial products 
of two vectors p, q. 

The terms along the vector r in these expressions do not deform the incom- 
pressible Earth and may be disregarded. Then ignoring the second-order terms, 
we can set W = grad W, the potential W at the right part being given by the 
expressions 

w = l -W^X^^;^'^''^^ , „ (2) 

[ - 2(1", ujy ~ [vixi + V2X2JX3UJ, r < Re, 

with the notations R, Re for the radii of the Earth and its fluid core, respectively. 

Adding the spherically symmetric term |r^a;^ to the right part (that does 

not affect the distribution of the density within the incompressible Earth) and 



10 



denoting cos 5 = {'p,U)/uj, p = r/r, we can present the potential W in each of 
these two domains within the Earth as the foUowing combinations of the zonal 
and tesseral harmonic functions: 



- |w-r^P° (cos S'), Rr<r<R, 

— |w^r^P2°(cos5) — (wia;i +1122:2)3^3'^, R < Re 



W^< !,,2.2p0 



In accordance with the general theory of Love numbers, action of the per- 
turbing spherical harmonics deforms the Earth's interior and the resulting de- 
formations induce the additional potential dW given on the Earth's spherical 
surface (of the mean radius R) by the expression 

dW\R = -'^uj^R'P°icoaS)~k^R\vipi+V2P2)p3t^, (3) 

o 

in which fc2, ^2 are the standard static and dynamic Love numbers, p = r/r = 
(pi, P2-, Ps) is the unit vector to the probing point out of the Earth. In the outer 
space, the first term at the right part generates the additional potential dWr 



dWr = -ifc2a;2-5-p0(cos5), r>R (4) 

which is the result of the tides aroused by rotation of the Earth as a whole, 
while the second term (proportional to A:|) brings about the tidal potential dWy 
caused by the differential rotation of the fluid core. The potential dWy has the 
form of the second-order tesseral harmonics: 



R^ 

dWv = -k^^-{vipi+V2P2)p3UJ, r> R. (5) 

For our aims, the zonal part of the potential Q may be ignored as it gen- 
erates nothing but the permanent tidal component of J2. Thus retaining only 
the tesseral terms (and keeping the old notation dWr) we obtain 

dWr = -k2^-{u}ipi+uj2P2)p3^, r > R. (6) 

Both these potentials, dW,. and dW^, should be added to the geopotential 
V which may be written in the standard form 



Gm.^ fR^ ^ 



V = — Y.\-) ^fc (^°s '^)(^fe ^°s-?'^ + *fc s^^-?'^)' (^) 



r ^ — ' \ r 

where S, A are the latitude and longitude of the probing point, P^ are associated 
Legendre functions. In order to present potentials ([5|) and ([6]) in the similar form, 

11 



the combinations pipa, P2P1, also should be expressed in terms of the associated 
Legendrc functions: 



P1/O3 = cos(5sin(5cos A = -P2^(cos(5) cos A, 

P2P3 = cos(5sin5sin A = -P2^(cos(5) sin A, 

o 

and then we obtain 



dWy = -k^{^—\ I— j P2^(cos(5)(wiCosA + W2sinA). 
Defining the so-called secular Love number kg in the standard way 

Gmh 

the expression for dW^ may be re-written in the form 

Cm / R\^ / k'^\ 
dW^=-r_f-J f-2 J j2PJ(cos^)(vicosA-^W2sinA)/w. (8) 

Such a presentation of the potential dWy demonstrates that its contribution 
into geopotential ([7]) has the form of corrections dcj , ds2 to the coefficients C2 , 
Sj- These corrections may be written in the simple form: 



The analogous expression is valid for the corresponding corrections dc^ , dsj 
proportional to ^2: 



In this relation we only account for the components due to the forced vari- 
ations of oJi, UJ2 (the so-called Oppolzer's terms) omitting the free Chandler's 
oscillation which is of much longer period. Hereinafter, the components of the 
near-diurnal oscillations in C2, S2 proportional to ^2 will be referred to as the 
Oppolzer's tides. 

One can see that the variables vi , V2 and uji , UJ2 in the right parts of relations 
([9|) and ([TO)l indeed oscillate with near-diurnal periods. For instance, the angular 
velocities wi , W2 in expression (jlOp may be presented in terms of the sidereal 
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time s and the time derivatives of the slowly changing angles of precession 
and nutation 9 by the Eulcr's kinematic relations: 



LOi — 6 coss + (j)sm9sms, 

UJ2 = —6sms + (psm6coss. (11) 

In analogous way, the components vi , V2 of the angular velocity of the fluid 
core may be transformed from the Earth-fixed frame into its components ni, n2 
in the inertial frame: 



Vi = Til cos s + 712 sm s 

V2 = — ni sins + n2 cos s. (12) 

Slowly changing variables ni, 712 can be expressed in terms of the nutational 
coefficients (and of the precession rate) making use of the differential equations 
of rotation of the deformable Earth with the fluid core; see (Moritz and Muller, 
1987). 

If dissipation in the fluid core is taken into consideration, the variables vi , V2 
should be calculated for the delayed time t — r. With all necessary accuracy, the 
linearized relations vi{t — t) = vi{t) — tvi and V2{t — t) = vi{t) — tv2 may be 
applied. Moreover, to calculate the time-delayed variables vi, V2, one can take 
into account only the time delay in the sidereal time s replacing in expression ([9]) 
vi by V1 — V25 and V2 by V2 + vi5 [5 = tw is the tidal phase lag). The dissipation 
brings about out-of-phase terms in the trigonometric presentation of the tidal 
effects. In the lERS standards, they are given for a theoretical value of the 
phase lag which, unfortunately, has not yet been verified by satellite data. 

Perturbing potential ([2]) of the centrifugal accelerations brings about not 
only the near-diurnal tidal variations ([3]), ([TU]) of the tesseral coefficients c\, s\ 
of the geopotential, but also near-diurnal displacements of sites: dH , dN, dE in 
the radial, northern and eastern directions, respectively. In our notations, the 
displacements caused by the fluid core may be written in the form 



dH = — i? I — 2- J J2 sin2(5(ui cos A -I- 'i;2 sinA)/a;, 

dN = -i? f -2- J J2Cos25(wiCosA + ?;isinA)/a;, (13) 

dE = —R 1 -2- I J2 sin5(— wi sin A -I- V2 cos A)/w, 

where A, 5 are the longitude and latitude of the site, /if, I2 are the dynamical 
Love numbers. 

The tidal displacements due to the centrifugal potential of the Earth as a 
whole may be obtained replacing the dynamical Love numbers /ijj '2 by the 
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static Love numbers /12, h and the angular velocities vi, vi of the differential 
rotation of the fluid core by the absolute angular velocities wi, ixii. Broken into 
the trigonometric series of the fundamental arguments, the corrections ([S]), (J13p 
are interpreted in the lERS standards (McCarthy, Petit, 2004) as the effect of 
the frequency-dependence of the Love numbers ^2 (or rather fc2i), /12 and li. 
It will be shown that the explicit analytical expressions of these trigonometric 
series have the following simple form (the proof is given at the end of Appendix) : 



dH = -i?f -^ j J2sin2(5^(C^cosA + 5'^sinA), 

dN = -i?Y-2- j J2cos2(5^(acosA + S'^sinA), (14) 

dE = -i?(-^) J2sin(5^(-asinA + S'^cosA), 



n)-"'--m^it 



A--' ■? 1: 7 I (15) 



where 



for 1/ ^ 0, 



S. = Qtsm{s + U) + Q-sm{s-,U), (16) 

a = Q+ cos(s + /„) + g; cos(s - /,) (17) 



•So = Qq sin(s) 
Co = Qocos(s), 

s is the sidereal time, f^, is the nutation argument of the frequency v in the 
nutational series given by the expressions 

9 = J^Se.cosf,, (18) 

ip = ^Sip^sinf^, (19) 

69^,Sipi, are coefficients of the nutations in inclination and longitude, and the 
coefficients Q'^, Q^ , Qq are expressed in terms of S9^,5ipi, by the relations 



Q+ = (69,- sin 9oS^p,)q,, (20) 

Q- = {d9, + sm9oSi^,)q^,, (21) 

Qo = — ^sin0o(l-K), (22) 
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with the notations 

ciu = ^^—{l-n), v^Q, (23) 



(24) 



p, a, Vc being the precession constant, the ratio of the moment of inertia of the 
fluid core to that of the Earth as a whole and the Free Core Nutation frequency, 
respectively. 

The corresponding effects of the Oppolzer's tides have the same form and 
may be obtained replacing /i2j ^21 ^2 by ^2, W, fe and making use of the following 
expressions for the coefhcients q^, Qq'. 

V 

Qo = sin 00- 

UJ 

We prefer the close form ([9]) , (fT3|) of the near-diurnal tidal effects instead of 
that recommended by lERS standards because, firstly, it is much simpler, and, 
secondly, this formulation explicitly demonstrates that the variations of c\, S2 
and the site displacements depend on the same combinations of the variables 
Vi, V2 (or wi, 1x12) and differ only by constant factors. Another advantage of 
this formalism written in the close form by relations ([9]), (fTO|). (fT3|) is that it 
may be applied not only with the analytical theories like lAU 2006 but also 
with numerical theories of the precession-nutation motion like that described 
in (Krasinsky, 2006), (Krasinsky, Vasyliev, 2006). This approach guarantees 
the consistency of calculations, unlike the formalism of lERS standards which 
suggests four different trigonometric series for four tidal effects under consider- 
ation. It is not clear how the recommended form of the near-diurnal tides (in 
which the Oppolzer's tides are not separated from the tides proportional to the 
dynamical Love numbers) may be applied varying numerically the dynamical 
Love numbers in order to calculate the corresponding partial derivatives. 

It is easy to see that the numerical values of the diurnal tides presented in the 
form of trigonometric series (|14|) -(|24 |) practically coincide with those given in 
the lERS standards (McCarthy, Petit, 2004). As an example, in Tab.l we give 
the amplitudes R^ (in mm) of the near-diurnal radial displacements calculated 
with the above expressions as well as the corresponding numerical values taken 
from lERS standards, see Tab. 7.5a of lERS Technical Notes 32 (McCarthy, 
Petit, 200). These displacements depend on the dynamical Love number /12 by 
relations 

(Rt]^_jij^li(Qt 

where Q+, Q^J are given by expressions p^ -p4 )) . and the value of ft.2 bas been 
taken in such a way that the largest amplitude (of the tidal constituent Ki 
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Tabic 3: CoefScients dRi, of radial displacements (mm) 



Ry 




lERS 


Period 


I 


/' 


F 


D 


n 


-0.04 


-0.04 


-0.08 


days 
9.13 


1 





2 





2 


-0.06 


-0.04 


-0.10 


13.66 








2 





1 


-0.28 
0.04 


-0.22 
0.02 


-0.51 
0.06 


13.70 
27.60 



1 






2 







2 



-0.05 


-0.01 


-0.06 


121.75 





-1 


-2 


2 


-2 


-1.20 


-0.10 


-1.23 


182.62 








2 


-2 


2 


-0.23 


-0.01 


-0.22 


6798.38 














-1 


11.72 


0.28 


12.00 


0.00 

















1.69 


0.04 


1.73 


-6798.38 














1 


-0.70 


0.00 


-0.50 


-365.30 





-1 











-0.14 


0.00 


-0.11 


-182.62 








-2 


2 


-2 



for which 1^=0) in the lERS standard would be equal to the sum of the corre- 
sponding terms proportional to /12 and /i2- Tab.l presents the coefficients Ri, 
giving not the frequencies v but by the corresponding periods (providing also 
the integer coefficients in the linear combinations of the fundamental arguments 
I, V , F, D, Q). As the periods are given with their signs, the upper index in 
the coefficients i?+ or R~ has been omitted. One can see that in the domain 
of the positive high frequencies ly, the contribution of the Oppolzcr's tides (the 
column rI ) is comparable with the effects of the fluid core. Tab.l demon- 
strates that the lERS coefficients indeed are the sums of the two components: 
those proportional to /cj and to fc2, the first component being the largest. (The 
coefficients of lERS and ours in the two last lines slightly differ). Analogous 
comparison with lERS standards were carried out also for the corrections dc\, 
ds\. The corresponding Fourier coefficients differ only by a constant factor 
from those of Tab. 4 and fixing properly this constant, a good agreement has 
been reached again. However, now we are interested in the absolute values of 
the coefficients which depend on the corresponding dynamical Love numbers. 
Serious discrepancies were found in this case. Indeed, for the maximal coeffi- 
cient (the constituent Ki, the zero frequency v) calculated with the a priory /cj 
we have obtained the value 589 x 10^^^ which considerably exceeds the value 
470 X 10-12 of jERS 2003. Making use of the new estimate k^ = 0.0595 de- 
rived in this work from the SLR data gives us the value 546 x 10" ^^ of this 
coefficient which diminishes the discrepancy but still significantly differs from 
the lERS value. These considerations reinforce our conclusion on the necessity 
of simultaneous estimating as J2 as fc2. 

Now we derive in brief the expressions given above for the coefficients Ci,, 
Si, of the near-diurnal tidal effects. We will express the pole coordinates wi,W2 
of the fluid core in terms of the coefficients of the adopted nutation 59, 6ip. For 
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the beginning, we restrict ourselves with the Poincare's model in which all tidal 
effects in wi, V2 are disregarded. As the adopted nutation theory has been con- 
structed in the framework of a more realistic Earth's rotation model which does 
account for the perturbations caused by the non-zero Love numbers, in this way 
the indirect perturbations of this type are automatically accounted for. Though 
a more refine approach might be developed in which the explicit dependence of 
Vi,i'2 on the Love numbers would be also considered, only the main effect of 
this sort is taken into account in (5+ because the resulting corrections are small 
enough. The following notations will be used: 

1. u — ui -\- iu2 and v — vi + iv2 are complex coordinates of the pole of the 
Earth and that of its fluid core, 

2. w is the angular velocity of the Earth, 

3. A, C are equatorial and polar moments of inertia. 

4. Ac,Cc are those of the core, 

5. e — {C — A)/ A, Be = (Cc ~ Ac)/ Ac are dynamic flatness of the Earth and 
that of its fluid core, 

6. a = Ac/ A. 

In the Poincare's approximation, the differential equations describing the 
time behavior of u, v may be written in the form (Poincare, 1987): 

u — iueuj — ail — iav + L (25) 

V = —iv{l + ec)uj — ii, (26) 

where i = %/— T and L = Li -\- iLi is a complex variable presenting the torques 
caused by outer forces. Transforming these equations into the normal form in 
which the time-derivatives are in the left side, we obtain 



u = lu Lo -\- IV w H , (27) 

1 — a 1 — a 1 — a 

f &c \ e L . , 

V — —iv 1 H (jj — iu Lj . (28) 

Here the torque L depends on the three Euler's angles: (the angle of 
nutation), -0 (the angle of precession) and (/) (the angle of the axial rotation). 
We can identify the variable (p with Greenwich Sidereal Time s. 

Supplementing these equations with the Euler's kinematic relations (jlip 
rewritten in the form 



ui = {ip sa\ 9 sn\ s -\- 9 cos s)/uj, 
U2 — (V'sin^coss — 0sins)/aj, 
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we obtain a close system of differential equations. The Euler's c relations may 
be presented in the equivalent complex form; 



u = Dexp{-is), (29) 

where 

D = {e + iijsme)/u}. (30) 

The torque L in equations (P7|) . (|28l) depends on the Euler's angles in the 
following way: 

L = cxp{-is)M{e,'ip) (31) 

where the function M is independent of the sidereal time s. 
Defining new variable Q by the relation 



exp(-is)g (32) 



we can consider D,Q as new variables and express differential equations (|27p. 
(|28|) in terms of these variables. Calculating the time derivatives of the both 
sides of identity (1^^ . we may ignore the time-dependence of the coefficient D. 
In the same approximation, s « w and we have 



ii = (— iwZ? + Z?) exp(— is) ~ — iwl?exp(— is). 

With definition (PTj) . equations (I^T]) . (^5)) may be written now in the following 
form: 



D = -Q_^+z_, 
1 — a uj 

Q = -iQ-^ iD- M. 

1 — a 1 — a 

Eliminating the complex function M from these equations, we obtain (ne- 
glecting the terms of the second order) the sought-for relation that ties the 
variables Q and D 



Q — —iQvc + iLoD, (33) 

with the notation Vc for Free Core Nutation frequency: 

Let 5N = 69 + iSijj sin ^o be the complex nutation variable. Due to identity 
((30|) . the variables SN and D are connected by the relation 
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D = — -iPsmeo, (34) 

UJ LO 

where the precession parameter p > is the coefficient of the linear trend in 
the precession angle ^ (we use the right-hand coordinate frame in which the 
precession motion is negative). The nutations 59, Sip may be presented in the 
form of the trigonometric series ^TE\\ . (IT5|) . In fact, the adopted nutations are 
referred to the instant coordinate frame, while we need the nutations in the fixed 
system of J2000. With the accuracy needed to calculate the small near-diurnal 
tidal effects, the difference may be ignored. The approximations f}^ = i' may be 
used due to the same reason. Now we present the complex nutational variable 
6N in the form of the trigonometric series: 



SN^Yl i^^ e^P(*/-) + ^-" exp(-*/.)) , (35) 



where 



N+ ^ ^{6e, + sm9o5iP), (36) 

N- = ^iS0,-smeoSyj). (37) 

The solution of equation (j33p for Q may be constructed in the similar form: 



= E (Q- exp(z/,) + Q- exp(-i/,)) + Qo- (38) 

Due to equations ([M|) - ([57)) . we can write the coefficients Q+, Qj7 and Qq as 



(39) 
(40) 

(41) 

Substituting relation (I5n])-(|1T|) into dTU]), (US]), we obtain the tidal effects 
under study in the form of relations P^ - (IM1) for the Poincare's model (for 
which the constant k defined by relation (IMl) is equal to zero). The Poincare's 
model is the first approximation to the adopted model based on the so-called 
SOS equations by Sasao, Okubo and Saito; see (Moritz and Muller, 1987). 
The next approximation (which is sufficient for our aims) may be obtained by 



Qt = 


ZZ + J/c 


5; = 


- Z ^^"' 


Qo = 


P 

— sin 6q — . 



19 



adding to the right part of equation ([^ the term kL. Repeating the above 
transformations with the modified in this way equations (|25l) . ((26|) . we obtain 
the analytical expressions (l^n|) - (P^ given earlier without proof. In fact, for 
calculations of the present paper we have used not these analytical expressions 
but more accurate values ni, n2 of the angular velocity of the fluid core (in the 
inertial coordinate frame) obtained by numerical integration of the SOS-type 
differential equations described in (Krasinsky, 2006). The analytical expressions 
were used just for a control, in order to be sure that the calculated corrections 
are close to those given in lERS standards. 
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